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ABSTRACT 

Diffuse interstellar bands (DIBs) trace warm neutral and weakly-ionized diffuse inter¬ 
stellar medium (ISM). Here we present a dedicated, high signal-to-noise spectroscopic 
survey of two of the strongest DIBs, at 5780 and 5797 A, in optical spectra of 666 
early-type stars in the Small and Large Magellanic Clouds, along with measurements 
of the atomic NaiD and CaiiK lines. The resulting maps show for the first time the 
distribution of DIB carriers across large swathes of galaxies, as well as the foreground 
Milky Way ISM. We confirm the association of the 5797 A DIB with neutral gas, and 
the 5780 A DIB with more translucent gas, generally tracing the star-forming regions 
within the Magellanic Clouds. Likewise, the NaiD line traces the denser ISM whereas 
the CaiiK line traces the more diffuse, warmer gas. The CaiiK line has an additional 
component at ^ 200-220 km s“^ seen towards both Magellanic Clouds; this may be 
associated with a pan-Magellanic halo. Both the atomic lines and DIBs show sub-pc- 
scale structure in the Galactic foreground absorption; the 5780 and 5797 A DIBs show 
very little correlation on these small scales, as do the CaiiK and NaiD lines. This 
suggests that good correlations between the 5780 and 5797 A DIBs, or between CaiiK 
and NaiD, arise from the superposition of multiple interstellar structures. Similarity 
in behaviour between DIBs and Nai in the SMC, LMC and Milky Way suggests the 
abundance of DIB carriers scales in proportion to metallicity. 

Key words: galaxies: ISM - ISM: atoms - ISM: lines and bands - ISM: molecules - 
ISM: structure - Magellanic Clouds 


1 INTRODUCTION 

The interstellar medium (ISM) plays a critical role in the 
star formation and chemical evolution of galaxies. Yet a full 
understanding remains elusive, with it being seen only in 
part or under special conditions. This is true in particular 
for the molecular component, as molecular hydrogen is virtu¬ 
ally invisible, challenging our ability to investigate how and 
when atomic gas becomes molecular. Even though the exact 
nature of their presumed molecular carriers is unknowrQ, 
the diffuse interstellar bands (DIBs) they generally cause in 
the spectra of stars are a powerful probe of the molecular 
ISM and the atomic-molecular interfaces. 

The Magellanic Clouds are two nearby (« 60 and 50 


^ Campbell et al. (2015) have recently reported confirmation of 
two near-IR diffuse bands as due to gas-phase C^q. 


kpc for the SMC and LMC, respectively) interacting, gas- 
rich dwarf galaxies with sub-solar ISM metallicities {Z ~ ^ 
and i Zq for the SMC and LMC, respectively). The reces¬ 
sion velocities (~ 150 km s“^ for the SMC and 250-300 km 
s“^ for the LMC) result in Doppler shifts with respect to ab¬ 
sorption in the Milky Way of ~ 3.5 A in the SMC and ~ 5.5 
A in the LMC at visible wavelengths. While some DIBs are 
broader than this, many are narrow enough to separate the 
Magellanic from the foreground Galactic absorption (Jen- 
niskens & Desert 1994). This is the case particularly for the 
5780 and 5797 A DIBs, which are intrinsically strong and 
their ratio is a clear diagnostic of exposure to UV irradia¬ 
tion and/or high electron density. 

Hutchings (1966) was the first to report a DIB in the 
Magellanic Clouds, viz. the broad 4430 A band. Subse¬ 
quent detections and upper limits were obtained by Blades 
& Madore (1979) and Houziaux, Nandy & Morgan (1980, 
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1985). Ehrenfreund et al. (2002) first detected the 5780 and 
5797 A DIBs in the SMC. While the UV radiation field seems 
to be the dominant factor affecting the relative strengths of 
DIBs in the Magellanic Clouds (Cox et al. 2006), other local 
conditions matter too (Welty et al. 2006), with Cox et al. 
(2007) highlighting metallicity effects both through reduced 
shielding against UV light and possibly directly through ele¬ 
mental depletions. Following pioneering efforts by van Loon 
et al. (2009) to map DIBs across the sky, van Loon et al. 
(2013) presented such maps towards the Tarantula Nebula 
in the LMC, showing a clear relation of DIBs with the diffuse 
ISM and with the radiation field. 

Here we present high signal-to-noise spectra from a ded¬ 
icated campaign to map the 5780 and 5797 A DIBs and the 
CallK and NalD lines across most of the SMC and a large 
part of the LMC. These are the first maps of this kind, 
revealing how the DIB carriers relate to the neutral and 
weakly ionized gas, the diffuse and star-forming ISM, and 
UV-shielded and irradiated conditions in these metal-poor 
environments. At the same time, we also map the Galac¬ 
tic foreground absorption, revealing small-scale structure in 
all these tracers, as well as CallK absorption arising in 
intermediate- and high-velocity clouds which are probably 
located in the Galactic Halo. 


2 DATA 

2.1 Observations 

The observations were carried out with the 3.9-m Anglo- 
Australian Telescope (AAT) at Siding Spring, Australia, 
during 3-6 November 2011 (programmes ATAC/lOB/51 and 
AAT/llB/17; PI: J.Th. van Loon). The two-degree field 
(2dF) 400-fibre positioning system was used in conjunc¬ 
tion with the AAH spectrograph mounted at the Coude 
focus, benefitting from a wide-held corrector and an at¬ 
mospheric dispersion compensator. The spectrograph has a 
dichroic beam-splitter allowing simultaneous observation of 
a “blue” and a “red” spectrum created with volume phase 
holographic (VPH) rehection gratings. Each spectrum was 
sampled with an e2v CCD with 2k pixels along the spectral 
and 4k pixels along the spatial direction. 

We used the 3200B and 2000R gratings with the 570-nm 
dichroic, with the central wavelength set to 3935 A for the 
blue arm and 5840 A for the red arm. This yielded a spectral 
resolving power of i? = A/AA ~ 8, 000 (38 km s“^), covering 
wavelength ranges of 3779-4084 A in the central hbre and 
3751-4054 A in the end fibres in the blue, and 5569-6101 A 
in the central fibre and 5528-6056 A in the end fibres in the 
red. 

The SMC field was observed in the first half of the night 
and the LMC field in the second. The weather conditions 
deteriorated throughout the observing run, with seeing from 
l!'3 to ~ 3” and increasing levels of cloud cover. Eventually, 
the SMC and LMC fields received total integration times 
of 12.5 hr and 13 hr, respectively, all split into individual 
30-minute integrations. 

Before and after each night the following calibration 
measurements were performed: bias (electronic offset), dark 
current level, arc lamp exposures for wavelength calibration, 
and fiat field observations on the sky and dome - the latter 


employing 20-W lamps for the red and 75-W lamps for the 
blue. 


2.2 Target selection 

The targets, fiducial stars (for acquisition) and sky positions 
are shown in figure 1; they were selected from the Magellanic 
Clouds Photometric Survey (Zaritsky et al. 2002). While the 
best targets for the study of interstellar absorption are early- 
type stars, which are bright and have relatively few and weak 
spectral features, their colours may have been reddened by 
interstellar dust and we must avoid too heavy a bias against 
such sight-lines. Therefore, initially targets were selected to 
have 13.5 < U < 15 mag, and —1 < (B — V) < 0.5 mag 
for the SMC but —1 < (B — V) < 0.4 mag for the LMC. 
A higher priority was given to stars with a reddening-free 
colour (Johnson 1963; Massey 2007) Q = (U — B) — 0.72 x 
{B — V) < —0.4 mag (red symbols in Fig. 1). When plotted 
using also the reddening-free Wesenheit index W — U —3.1 x 
[B — V) (Madore 1982, using the canonical value Rv = 3.1) 
the early- and later-type stars separate much more clearly 
than in the unreddened colour-magnitude diagram. Indeed, 
the early-type stars trace the star forming regions whereas 
the later-type stars are concentrated towards the LMC bar 
in the South-West corner of the 2° field (Fig. 1). The later- 
type stars, while not ideal targets for our study, were not 
excluded from selection entirely as they could still fill some 
fibres that could not otherwise be allocated. In Appendix 
A we describe the distribution over known spectral types, 
confirming the success of our selection strategy. Fiducials 
(green triangles in Fig. 1) were selected from 13th-mag (V) 
stars with 0.3 < {B — V) < 0.6 mag. Targets or fiducial 
stars were rejected if other stars were known to affect the 
signal within the fibre significantly. Sky positions (blue dots 
in Fig. 1) were chosen to be clear of known stars within or 
near to the fibre. 

The CONFIGURE software was used to allocate fibres to 
targets, up to 8 fiducial stars, and 25 empty positions for 
sky measurements to be subtracted from the target stars’ 
spectra. An annealing scheme was employed to optimise the 
target allocation, yielding « 350 targets per 2dF field. The 
final selection of observed targets - 338 in the SMC and 
328 in the LMC - is shown in figure 2 and table 1 (made 
available in full at CDS). 

2.3 Data processing 

The data were processed using the 2dFdr tool (Lewis et al. 
2002). The reduction steps include: bias subtraction, dark 
subtraction, fiat-fielding, tram-line mapping of the (curved) 
spectra on the CCD, spectra extraction, arc lamp spectral 
line identification, wavelength calibration, fibre throughput 
calibration and sky subtraction. The blue and red spectra 
are processed separately. 

Accurate subtraction of sky emission lines can be a chal¬ 
lenge for fibre-fed spectrographs. While the standard sky 
subtraction in 2dFdr produced good results, the very strong 
telluric NalD lines were not always removed entirely (see 
also the high resolution fibre-fed spectra in van Loon et al. 
2013). The residual telluric emission was accounted for at 
the analysis stage (see next section). 
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Figure 1. Distribution of SMC {top row) and LMC {bottom row) targets (high priority: red; low priority; magenta), fiducials (green 
triangles) and sky positions (blue dots) on the sky {left), in optical colour and brightness {middle) and in reddening free indices {right). 


Table 1. List of observed targets. Only the first three objects in each of the SMC and LMC are shown here. The ID numbers are from 
our own target lists. Both equatorial and Galactic coordinates are given as the latter are useful in interpreting the foreground absorption. 
Spectral types and luminosity classes are based on a literature search and only available for 174 SMC and 61 LMC targets (see Appendix 
A). B- and V-band photometry are from Zaritsky et al. (2002). The full version of this table is available on CDS. 


ID 

RA (*' “ =) 

Dec (° ' ") 

1 (°) 

6(°) 

B (mag) 

V (mag) spectral type 

SMC: 

983 

1 05 49.02 

-72 48 18.0 

301.44836 

-44.28177 

14.58 

14.79 

942 

1 04 05.97 

-72 43 50.9 

301.61804 

-44.36510 

14.02 

14.16 

762 

0 58 11.26 

-72 36 10.0 

302.22449 

-44.51616 

14.36 

14.21 

LMC: 

1856 

5 41 16.56 

-69 14 02.2 

279.59109 

-31.43029 

14.31 

14.20 

1913 

5 41 59.34 

-69 14 56.5 

279.60104 

-31.36586 

15.11 

14.97 

1281 

5 36 09.39 

-69 13 58.0 

279.65109 

-31.88117 

14.58 

14.50 


The spectra from all of the nights were processed indi¬ 
vidually but then combined using custom-written routines 
in IDL: first they were rectified and normalised, then the 
root-mean-square noise level (standard deviation in a line- 
free region) was determined, the reciprocal value of which 
acted as the weight in averaging the spectra. 


3 ANALYSIS 

In this work we concentrate on the ISM tracers in the spec¬ 
tra, viz. the CallK (the H component 34.8 A further to the 


red is obliterated by He), NalDi and D 2 , and the 5780 A 
and 5797 A DIBs, both for the SMC and LMC as well as 
the Galactic foreground. 

The spectral regions around the above absorption fea¬ 
tures were fitted using a combination of a first-order poly¬ 
nomial (baseline) and Gaussian functions. The actual fitting 
was coded in idl using the mpfit routines. 

The constraints on the fitting are summarised in table 2. 
In particular, the wavelength of the 5797 A DIB was guided 
by that of the stronger 5780 A DIB and allowed to vary 
less to avoid spurious results. Likewise, the wavelength and 
width of the Di (red) component of the Nal doublet were 
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Figure 2. Distribution of observed targets in the SMC (top) and 
LMC (bottom). Some features of interest are indicated, including 
30 Doradus, SN 1987A, etc. 


tied to that of the stronger Dg (blue) component - except in 
the LMC where the opposite strategy was used because its 
D 2 line coincides with the Galactic foreground Di line. Over¬ 
lapping components were fitted simultaneously. The peaks 
of the Gaussian functions were forced to be negative, but an 
additional positive-valued double-Gaussian function (with a 
fixed separation of 5.97 A) was included around the rest 
wavelengths of the Nal D lines to account for telluric emis¬ 
sion. It was noted that the SMG and LMC components of 
the Ca II K line often showed a second component. Towards 


Table 2. Overview of the Gaussian fitting constraints, in A. 


Feature 





^0 

guess 

range 

guess 

range 

SMC, internal: 





CallK (blue) 

3935.6 

± 2.00 

0.2 

0.10 -0.40 

CallK (red) 

3936.6 

± 1.00 

0.2 

0.10 -0.40 

5780 DIB 

5783.6 

± 1.00 

0.6 

0.036-0.90 

5797 DIB 

A578O+I6.5 

± 0.75 

0.5 

0.25 -1.00 

Nal D 2 

5893.0 

± 2.00 

0.5 

0.25 -1.00 

Nai Di 

Ad2 + 5.97 


<^02 


Galactic, in SMC direction: 




CaiiK 

3933.6 

± 2.00 

0.3 

0.15 -0.60 

5780 DIB 

5780.1 

± 1.00 

0.6 

0.036-0.90 

5797 DIB 

A578O+I6.5 

± 0.75 

0.5 

0.25 -1.00 

Nal D 2 

5890.0 

± 2.00 

0.5 

0.25 -1.00 

Nal Di 

Ad2 + 5.97 


( 7 D 2 


LMC, internal: 





CallK (blue) 

3937.0 

± 1.00 

0.2 

0.12 -0.30 

Ca II K (red) 

3937.5 

±0.50 

0.2 

0.12 -0.30 

5780 DIB 

5786.0 

±1.00 

0.6 

0.036-0.90 

5797 DIB 

A578O+I6.5 

±0.75 

0.5 

0.25 -1.00 

Nal D 2 

Adi - 5.97 




Nal Di 

5901.0 

±2.00 

0.5 

0.25 -1.00 

Galactic, in LMC direction: 




Ca IIK (blue) 

3933.6 

±1.00 

0.2 

0.12 -0.30 

Ca II K (red) 

3935.0 

±0.50 

0.2 

0.12 -0.30 

5780 DIB 

5780.8 

±1.00 

0.6 

0.036-0.90 

5797 DIB 

A578O+I6.5 

±0.75 

0.5 

0.25 -1.00 

Nal D 2 

5890.0 

±2.00 

0.5 

0.25 -1.00 

Nal Di 

Ad2 + 5.97 


( 7 D 2 



the LMC, a second component was also seen in the Galactic 
foreground. 

The equivalent width EW was determined by integrat¬ 
ing the Gaussian function: 

EW = J A exp 1^- = AaV^, (1) 

where A is the peak intensity, Aq the central wavelength 
and a the width (with the Full Width at Half Maximum, 
FWHM = 2.355(7). The errors eg resulting from mpfit- 
EXPR were scaled according to the deviations from the 
ht and the degrees of freedom / (the number of spectral 
elements minus the number of htting parameters) as in: 

e = eos/x^/f (2) 

Hence the error in the equivalent width is obtained from: 

e[EW] = \J2'K ((e[A]a)^ + (Ae[(7])2) (3) 

All spectra and fits were vetted by eye, guided by the 
value for FTT/e[i5VK], and only fits deemed reliable were re¬ 
tained. Rejected data include a few spectra with broadened 
CallK absorption arising in the stellar photosphere (Ap¬ 
pendix A presents an overview of the spectral types and their 
mostly negligible effect on the measured equivalent widths). 
The full tables of measurements are available at GDS, with 
the first three entries shown in tables 3 and 4. 

Three of our targets (all in the LMC) have also been 
subject of measurement of the 5780 A DIB in previous stud¬ 
ies. Our target 1491 (BI253) was measured in a high res¬ 
olution spectrum by Welty et al. (2006); their EW 57 S 0 = 
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Table 3. Spectral line measurements for SMC sight-lines. EW 
stands for equivalent width; e[EW] its uncertainty. Flags / are 
1 or 0 for a detection or non-detection, respectively; The full, 
transposed version of this table is available on CDS. 


ID 

983 

942 

762 

SMC, internal: 




EIFcallKblue (A) 

0.224 

0.185 

0.310 

e[EVF] Call K blue (A) 

0.019 

0.005 

0.006 

(ilIF/e[i?IV])caIIKblue 

11.6 

36.1 

48.5 

EWcallKred (A) 

0.040 

0.041 

0.154 

e[EVF] Call K red (A) 

0.012 

0.004 

0.005 

(EIF/e[i?lF])caIIKred 

3.4 

9.1 

30.5 

EW 578 O (A) 

- 

0.006 

0.052 

e[EW]s,7so (A) 

- 

0.001 

0.002 

{EW/e[EW])5780 

- 

6.5 

23.3 

EW 5797 (A) 

0.016 

- 

0.009 

e[EVF]5797 (A) 

0.007 

- 

0.001 

iEW/e[EW])s797 

2.3 

- 

10.7 

ElLNaID2 (A) 

0.096 

0.096 

0.246 

EWNalDl (A) 

0.129 

0.064 

0.164 

e[EVF]NalD2 (A) 

0.008 

0.004 

0.005 

e[EVF]NaIDl (A) 

0.008 

0.004 

0.004 

(EIF/e[ii;iT])NaID2 

12.2 

21.7 

46.2 

(ElT/e[EVF])NaiDi 

16.0 

18.0 

38.2 

fCall K blue 

1 

1 

1 

/5780 

0 

0 

1 

/5797 

0 

0 

1 

/NalD 

0 

1 

1 

Galactic, in SMC direction: 



EVFcallK (A) 

0.190 

0.126 

0.088 

e[EVF]caiiK (A) 

0.019 

0.004 

0.005 

(EIF/e[i?lT])caiiK 

10.0 

31.4 

19.0 

EIT 5780 (A) 

0.033 

0.024 

0.022 

e[EW]57so (A) 

0.005 

0.002 

0.003 

iEW/elEW])s780 

6.7 

13.7 

7.6 

EIT 5797 (A) 

0.016 

0.013 

0.025 

e[EVF]5797 (A) 

0.002 

0.002 

0.002 

iEW/e[EW])s797 

6.8 

5.6 

10.9 

ElLNaID2 (A) 

0.136 

0.086 

0.209 

EIFNalDl (A) 

0.048 

0.056 

0.148 

e[EVF]NaID 2 (A) 

0.008 

0.006 

0.012 

e[EVF]NaIDl (A) 

0.005 

0.005 

0.011 

(ElT/e[ElT])NaID2 

16.1 

15.0 

17.8 

iEW/e[EW]),,^ioi 

9.2 

11.5 

13.9 

fCall K 

1 

1 

1 

/ 578 O 

0 

0 

0 

/S797 

0 

0 

1 

/nuID 

1 

1 

1 


0.075 ± 0.016 A agrees with our EW5780 = 0.056 ± 0.003 
A. Our targets 1611 (VFTS 532) and 1625 (VFTS604) had 
been measured in similar high resolution spectra by van 
Loon et al. (2013); they got EW 578 O = 0.175 ± 0.007 A 
for star 1611 compared to our EW 57 S 0 = 0.194 ± 0.006 A, 
and EW 878 Q = 0.123 ± 0.009 A for star 1625 compared 
to our EW 57 S 0 = 0.137 ± 0.003 A. All three comparisons 
are consistent within 2-3cr. Reassuringly, the measurements 
of the sodium lines also agreed very well: for star 1611, 
7?lFNaiDi = 0.463 and 0.466 A in van Loon et al. (2013) 
and this work, respectively, while for star 1625 these values 
were FiVFNaiDi = 0.467 and 0.451 A, respectively. 


Table 4. As table 3 but for LMC sight-lines. 


ID 

1856 

1913 

1281 

LMC, internal: 




EITcallKblue (A) 

0.243 

0.255 

0.035 

e[EVF] Call K blue (A) 

0.005 

0.007 

0.003 

(EW/e[EW])c^UKhiu. 

45.1 

37.3 

11.3 

ElFcallKred (A) 

0.169 

0.190 

0.168 

e[EVF] Call K red (A) 

0.005 

0.009 

0.005 

(EIF/e[i?lF])caIIKred 

31.4 

21.2 

32.2 

EW 5780 (A) 

0.141 

0.015 

0.051 

e[EIF]5780 (A) 

0.003 

0.003 

0.003 

iEW/elEW])s780 

55.5 

4.8 

20.2 

EW5797 (A) 

0.030 

0.017 

0.012 

e[EVF]5797 (A) 

0.003 

0.002 

0.004 

iEW/e[EW])s797 

11.8 

10.2 

3.5 

EWNaID2 (A) 

0.544 

0.524 

0.308 

EWNalDl (A) 

0.424 

0.216 

0.232 

e[EVF]NaID2 (A) 

0.013 

0.049 

0.017 

e[EVF]NaIDl (A) 

0.006 

0.005 

0.005 

(EIF/e[ii;VF])NaID2 

41.6 

10.8 

17.6 

(EIF/e[EVF])NaiDi 

66.5 

39.7 

46.1 

fCall K blue 

1 

1 

1 

/5780 

1 

0 

1 

/5797 

1 

0 

0 

/NalD 

1 

1 

1 

Galactic, in LMC direction: 



ERCallKblue (A) 

0.128 

0.118 

0.127 

e[ElF] Call K blue (A) 

0.005 

0.004 

0.006 

(ilIF/e[EIV])caIIKblue 

28.5 

30.6 

22.6 

ElFcallKred (A) 

0.034 

0.026 

0.121 

e[ElF] Call K red (A) 

0.005 

0.006 

0.004 

(£;iF/e[EIT])caIIKred 

7.5 

4.6 

33.7 

EW 578 O (A) 

0.017 

0.027 

0.010 

e[EW]8780 (A) 

0.001 

0.004 

0.001 

{EW/e[EW])8780 

11.6 

6.6 

9.5 

EWs797 (A) 

0.009 

0.016 

0.005 

e[EVF]5797 (A) 

0.001 

0.004 

0.002 

{EW/e[EW])5797 

7.2 

4.3 

2.9 

EIFNaID2 (A) 

0.195 

0.171 

0.312 

WNalDl (A) 

0.254 

0.099 

0.257 

e[EVF]NaID2 (A) 

0.041 

0.019 

0.186 

e[EVF]NaIDl (A) 

0.061 

0.015 

0.154 

(ElF/e[ii;VF])NaID2 

4.7 

9.0 

1.7 

(EVF/e[EVF])NaiDi 

4.2 

6.4 

1.7 

fCall K blue 

1 

1 

1 

/5780 

1 

1 

1 

fh797 

0 

0 

0 

/NalD 

1 

1 

1 


4 RESULTS 

The average line profiles are shown in figure 3 (with the fits 
overplotted in red). From the fitted line shape, two veloc¬ 
ity components of the CallK line are clearly visible in the 
Galactic foreground towards the LMC (Fig. 3e). Many spec¬ 
tra recorded towards individual stars also showed a second 
component internal to the SMC and LMC. The redshift of 
the LMC is large enough for its 5780 A DIB to be separated 
from that of the Galactic foreground (Fig. 3f). Towards the 
SMC the two velocity components are blended (Fig. 3b), 
but simultaneous fitting of two Gaussian functions worked 
well in accounting for their respective contributions. The 
narrower 5797 A DIB is resolved both towards the LMC 
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Figure 3. Average spectra for the SMC {top row) and LMC {bottom row) around the {from left to right:) CallK, 5780 A DIB, 5797 A 
DIB and NalD lines. The average fits are overplotted in red. 


(Fig. 3g) and towards the SMC (Fig. 3c). The NaiD lines 
are well separated for the SMC (Fig. 3d) but towards the 
LMC the LMC D 2 component is close to the Di component 
of the Galactic foregronnd absorption (Fig. 3h). But again, 
simultaneons fitting worked well in resolving the two compo¬ 
nents. This is confirmed in the maps we show later, which - 
as expected - are very different for the different components. 

While the spectral resolving power - and intrinsic width 
of the DIBs - did not allow us to identify multiple kine¬ 
matic components other than (sometimes) a second one in 
the CallK, it is very likely that multiple kinematic compo¬ 
nents are present in all absorption features and sight-lines 
(cf. van Loon et al. 2013). Likewise, the CallK and NalD 
absorption may be saturated in some of these components, 
and an increase of their equivalent width can be caused by 
stronger line wings as well as additional kinematic compo¬ 
nents. The Di line has half the intrinsic strength of that of 
the D 2 line and shonld therefore be less affected by satura¬ 
tion effects. While we caution against the use of the equiv¬ 
alent widths of the CallK and NalD lines in an absolute 
sense, a comparative analysis can still be meaningful. 

Remarkably, the Ca ll K absorption towards the SMC is 
relatively strong compared to the Galactic foreground and 
LMC sight-lines; the DIBs are considerably weaker, though 
it must be said that the bulk of the absorption within the 
SMC is concentrated in areas much smaller than the 2° field 
and individual sight-lines can show quite strong DIBs. 


4.1 Correlations 

Figures 4-6 show plots of the EW of the 5780 vs. 5797 A 
DIBs, 5780 A DIB vs. NaiDi and 5797 A DIB vs. NaiDi, 
respectively. (Almost) no correlation was discernible be¬ 
tween the DIBs and CallK. Plotted are not only the secure 
detections in both relevant features, but also (with different 
colours/symbols) less certain detections and non-detections. 
While some outliers and some of the scatter are due to un¬ 
certainties, there are clear positive correlations between the 
5780 and 5797 A DIBs and with NalD but with real (in¬ 
trinsic) scatter also. The equivalent width of the NalDi line 
grows slower than in proportion to column density because 


of saturation effects, giving the false impression of acceler¬ 
ated growth of the DIBs. 

The Galactic foreground shows little difference between 
the SMG and LMC directions, but the SMC itself clearly 
has far fewer sight-lines with strong DIB absorption than 
the LMC (Fig. 4). Strong DIB absorption is seen internal 
to the SMC towards star 135 (SMC 10196 in Massey 2002) 
in the centre of the South-Western star-forming complex: 
F/IT 5797 = 0.079 A and EW^rso ~ 0.280 A. The NaiD 
absorption is also strong in that direction. This sight-line 
passes in between the LHA-115N25 OB assocation and the 
SMC-B2 3 molecular cloud; possibly it crosses the irradi¬ 
ated skin of a dense cloud, which may explain the strong 
DIB absorption. Similar strength of DIB absorption is seen 
at several locations within the LMC; the correlation appears 
to be saturating at EWsno ~ 0.30 A while the 5797 A DIB 
continues to strengthen beyond EW 5797 > 0.1 A. The two 
stars that show the strongest 5797 A DIB in figure 4 are 
both located within the Tarantula Nebula (cf. van Loon et 
al. 2013). 

In Table 5 the results are summarised from a linear 
regression analysis and a Student t test for the 5780 and 
5797 A DIBs: the correlation is significant if t > fcriticai, 
where 


t = r 


N-2 

1 _ J.2 


( 4 ) 


with correlation coefficient r and N data points {N — 2 de¬ 
grees of freedom in this case). Hence we conclude that a good 
DIB-DIB correlation is seen within the SMC and LMC but 
that the Galactic foreground is probably too transparent to 
show a significant correlation. The 5780/5797 DIB ratio is 
very similar between the SMC, LMC and the Local Bubble 
- where EW5780 = 0.023(±0.005) -f 2.48(±0.20) x EW5797 
(Bailey et al. 2015). There is a hint that the relation 
is steeper (relatively stronger 5780 A DIB) in the SMC, 
which could be explained as due to a harsher radiation field 
from metal-poor stars (reduced line-blanketing) in a metal- 
poor ISM (reduced attenuation by dust), but this is only 
marginally significant. 

Both the DIBs and NalD are weaker in the SMC 
than in the LMC (Figs. 5 & 6 ), but they behave simi- 
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Figure 4. The equivalent width of the 5780 and 5797 A DIBs for the SMC {top left), LMC {top right) and the Galactic foreground in 
the direction of the SMC {bottom left) and LMC {bottom right). 


Table 5 . Linear regression analysis of DIB-DIB correlations of the form -BW5780 = a + b x -EIT5797 with correlation coefficient r. Where 
the Student t statistic is given, the critical value of t is for a 1% chance that the null hypothesis is true. 


location 

N 

a (A) 

b 

r 

t 

^critical 

verdict 

SMC, internal 

16 

0.007 ± 0.005 

2.93 ± 0.47 

0.86 

6.3 

2.6 

significant 

LMC, internal 

44 

0.036 ± 0.013 

2.17 ±0.28 

0.77 

7.8 

2.4 

significant 

Galactic, in SMC direction 

38 

0.032 ± 0.005 

0.16 ± 0.31 

0.09 

0.5 

2.4 

insignificant 

Galactic, in LMG direction 

25 

0.016 ± 0.003 

1.24 ±0.53 

0.44 

2.3 

2.5 

very marginal 


larly. Both show what appears to be a threshold, or a more 
slowly picking up of DIB strength as Nal D increases, around 
i?WNaiDi ~ 0.1 A. A similar behaviour was noticed in the 
Local Bubble (Bailey et al. 2015); it probably is due at least 
in part to the NalD saturating. The Galactic ISM is clearly 
more transparent in front of the SMC than in front of the 
LMC, which is probably related to the larger distance to the 
Galactic plane of the former. Quite strong absorption is seen 
in sight-lines close to the LMC in the Local Bubble survey 
(Bailey et al. 2015) though no sight-lines were observed by 
these authors within 10° of the SMC. 


4.2 Maps 

We constructed maps of the 5780 A DIB, 5797 A DIB, 
CallK and NalD 2 equivalent widths within the SMC (Fig. 
7), within the LMC (Fig. 8) and within the Galactic fore¬ 
ground (Figs. 9 & 10). We do not show the NalDi maps 
but they look identical to the D 2 maps. We overlaid the 
positions of the stellar probes; as absorption was not al¬ 
ways detected the maps are generally more sparsely sam¬ 
pled than the distribution of targets - which are already not 
distributed uniformly across the sky. The maps are filled 
in using triangulation, which preserves the true information 
(measurements) and does nothing but interpolate between 
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Figure 5. The equivalent width of the 5780 A DIB versus that of the NalDi, for the SMC {top left), LMC {top right) and the Galactic 
foreground in the direction of the SMC {bottom left) and LMC {bottom right). 


measurements without smoothing. The non-detections (ma¬ 
genta symbols) provide a meaningful addition to the maps 
as they outline areas of very low absorption. The SMC maps 
of the 5780 and 5797 A DIBs are based on just 35 and 19 
firm detections, respectively; for the LMC these are 91 and 
48, respectively. The CallK line was detected in most sight¬ 
lines in both galaxies and the NalD in about half of the 
sight-lines. 

In the SMC (Fig. 7) the DIBs are strongly associated 
with the South-Western part of the SMC “bar”, where 
star formation is widespread. The EWsrso/EWsrQr ratio ap¬ 
pears to peak slightly offset from that region, perhaps trac¬ 
ing harsher conditions in the periphery of the star-forming 
molecular clouds as compared to those within the clouds. 
The South-Western part of the bar is also where NalD 
peaks though it is detected more globally across the SMC. 
No DIBs were detected around the mini-starburst N 66 fur¬ 
ther North in the bar. The blue velocity component of 
CallK is detected across the SMC, tracing the generally 
hot, widely distributed gas typical of the metal-poor ISM in 
the SMC, but the red velocity component is concentrated 
towards the Wing, to the South-East of the bar. The red 
component might be situated behind the bulk of the stars. 


hence the lack of detected absorption in this gas. There is 
no resemblance between the CallK maps and those of the 
DIBs, but it can be noted that where the DIBs are strongest 
the CallK is relatively weak. This shows a stronger affilia¬ 
tion of DIBs with Nal-bright material than with Call-bright 
material. 

The LMC maps (Fig. 8) confirm this behaviour: The 
DIBs are strongest in molecular clouds within star-forming 
regions, but not too close to recently formed OB associa¬ 
tions such as the central cluster R136 (30 Doradus) within 
the Tarantula Nebula. This was also found in the DIB maps 
of the Tarantula Nebula by van Loon et al. (2013). The 
EWs 7 so IEWz 7 Q 7 ratio is high right at the location of R136, 
as well as towards the West where the ISM becomes more 
diffuse (less self-shielded). The NalD is concentrated in the 
same regions of great DIB strength, though more spread 
out. The CallK data are interesting: as in the SMC, the 
blue component is seen across the LMC also in areas with 
weak Nal D, though the densest parts tend to be situated in 
the East. The red component of CallK resembles the NalD 
map more closely; here the East-West gradient is particu¬ 
larly evident, with the strongest absorption along the East¬ 
ern, leading edge of the LMC. Again we see a fair correspon- 
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Figure 6. The equivalent width of the 5797 A DIB versus that of the NalDi, for the SMC (top left), LMC (top right) and the Galactic 
foreground in the direction of the SMC (bottom left) and LMC (bottom right). 


dence between DIBs and NalD but not with CallK except 
for a few locations of generally high column density - e.g., 
the blue component of CallK in the South of the Tarantula 
Nebula and the red component to the West of it. 

The maps of the Galactic foreground (Figs. 9 & 10) 
bear no resemblance to those of the SMC and LMC, which 
lends credibility to both the Galactic and Magellanic maps 
we created. In both the SMC and LMC sight-lines the DIBs 
show gradients over scales of a degree, with possibly smaller 
scale structure - similar to what was seen in the extra-planar 
gas in front of the globular cluster uj Centauri by van Loon 
et al. (2009). For instance, the 5797 A DIB (Fig. 9) traces 
a coherent presumably fairly neutral cloud lying across the 
South-Eastern portion of the SMC field, whereas the 5780 A 
DIB identifies material towards the South-Western corner of 
the LMC field. The EW^rsol ratio is high in places 
where the 5797 A DIB is very weak, not where the 5780 
A DIB is strongest. This suggests that the ratio is high in 
relatively diffuse gas under special conditions. 

The foreground CallK and especially the NalD (Fig. 
10 ) is structured with high contrast on much smaller scales, 
~ Of 1. This makes it hard to compare these maps with those 
of the DIBs (Fig. 9), though where the NalD is strong 


(Southern half of the SMC field, around Dec ~ —69?7 in 
the LMC field) the stronger DIBs are seen also. It is possi¬ 
ble that the DIBs are structured on similar small scales as 
the atomic absorption, but the noise and fewer detections of 
the much weaker DIBs prevents direct proof. 

5 DISCUSSION 

5.1 The carriers of the DIBs 

The 5780 A DIB is approximately 2-3 times as strong as the 
5797 A DIB throughout most of the SMC and LMC. This 
contrasts with the very high ratio of ElLsrso/EW 5797 ~ 6 
found towards R136, the ionizing cluster of the Tarantula 
Nebula (van Loon et al. 2013). The latter is an extreme 
environment and hence not typical of the ambient ISM 
within the Magellanic Clouds. The EWNaiDi/EWsrso ~ 3- 
4, on the other hand, is very similar between our global 
maps and the R136 sight-line. The Local Bubble and 
its immediate surroundings comprise a typical range of 
EIF 5780 /EW 5797 ~ 2-5 (Farhang et al. 2015), concentrated 
around EWsrsoIEWsrgr ~ 3 (Bailey et al. 2015), which en¬ 
compass the ratios found in the Magellanic Clouds. Very 
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Figure 7. Maps of the SMC in the equivalent width of (top row:) the 5780 and 5797 A DIBs and their ratio, and (bottom row:) the blue 
and red CallK components and NalD 2 . The locations of the target stars are indicated with the small white squares (detections) and 
magenta diamonds (non-detections). 


high values of EWsrsojEWsrgr ~ 5-9 are seen within the 
Local Bubble (Bailey et al. 2015); the highest reliable ra¬ 
tios within the Magellanic Clouds are EW^rso /SIA 5797 ~ 7. 
The ratio i51LNaiDi/i?W5780 ~ 2-3 found around the Local 
Bubble (Farhang et al. 2015) is again commensurate with 
what is seen across the Magellanic Clouds. 

We note that the DIBs are conspicuously weak in the 
Southern molecular ridge South from 30 Doradus (towards 
N 158) which is traced by NaiD and (to a lesser degree) 
CallK. This was also noted by van Loon et al. (2013) and 
may indicate the disappearance (or lack of excitation) of 
DIB carriers in dense, UV-shielded environments. 

The ionization potential of Na° is 5.1 eV, while that of 
Ca° and Ca+ are 6.1 and 11.9 eV, respectively. Comparing 
this to the ionization potential of H°, 13.6 eV, it is clear 
that in Hll regions neither Nal nor Call correspond to the 
dominant stage of ionization. However, in cooler, generally 
neutral clouds Call will diminish with respect to Nal, which 
is further compounded due to depletion of calcium atoms 
into grains. It is not yet clear whether DIB carriers survive 
when immersed in photon baths of > 6 eV, so mapping 
the DIBs alongside the Nal and Call may elucidate this 
point. Our maps show that the 5780 and 5797 A DIBs more 
closely trace Nal than Call, with the 5797 A DIB tracing 


cooler and/or more shielded structures than the 5780 A DIB 
does. This suggests that the carrier of the 5780 A DIB is 
removed by particles/photons with energies in excess of 6 
eV, while the carrier of the 5797 A DIB is already removed 
at energies in excess of 5 eV. As the 5780 A DIB stops 
growing in tandem with the 5797 A DIB for large column 
densities, the carrier of the 5780 A DIB probably needs to 
be sustained by particles/photons with energies of a few eV, 
making it likely to be a cation. 

Remarkably, the equivalent-width diagrams of either 
the 5780 or 5797 A DIB versus NaiD (Figs. 5 & 6) are 
indistinguishable between the SMC and LMC but for the 
fewer sight-lines with strong absorption in any of these trac¬ 
ers (an observation which is facilitated by our decision to 
plot all diagrams between SMC, LMC and Galactic fore¬ 
ground on exactly the same scale). Figure 11 shows this 
for the 5780 A DIB, by comparing all good (> IOct) detec¬ 
tions between the SMC, LMC and Galactic foreground in 
one graph. This implies that any metallicity dependence af¬ 
fects both sodium and DIB carriers in essentially the same 
manner. This is very curious, especially if the carriers of the 
DIBs are complex molecules. In that case one would expect 
a stronger than linear relation between overall metal content 
and DIB carrier. As the sodium atomic abundance scales ap- 
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Figure 8. As figure 7, but for the LMC. 


proximately in proportion to metallicity, apparently so does 
the DIB carrier. It means, in turn, that if the DIB carriers 
are indeed multi-atomic compounds their formation must be 
efficient and largely complete and this material must be re¬ 
silient. For instance if they are carbon-based, for instance in 
the form of fullerenes, then much of the interstellar carbon 
must be locked within the DIB carriers. 


5.2 Structure of the Galactic and Magellanic ISM 

The Galactic foreground maps of 5780 and 5797 A DIBs 
and CallK and NalD absorption reveal structures on mul¬ 
tiple scales - up to a degree down to as small as ~ 0?1, 
which for typical distances d ~ 200 pc means linear scales 
of ~ 0.35-3 pc. The same has been seen in both DIBs 
and NalD through extra-planar sight-lines (van Loon et al. 
2009), and in CallK and NalD through intermediate- and 
high-velocity clouds (IVCs and HVCs, respectively; Smoker, 
Fox & Keenan 2015). Curiously, the correlation between the 
two DIBs, or Call and Nal, is very weak at best - the maps 
show no conspicuous relationship between each of these pairs 
of tracers. This suggests that where correlations do exist it 
is likely due to the cumulative effect of a superposition of 
structures in a multi-phase ISM in which each individual 
structure is dominated by one or the other. 


The leading edge of the LMC is dense in both Na° and 
Ca"*". The DIBs are strongest in the star forming regions to 
the South and West of the Tarantula Nebula; they clearly 
trace the warmer - but not strongly-irradiated - cloud lay¬ 
ers. This is different from both dense, UV-shielded molec¬ 
ular cloud tracers such as CO and far-IR emission on the 
one hand (or reddening - cf. Tatton et al. 2013), and photo¬ 
dissociation regions on the other hand. It would be inter¬ 
esting to see how well DIBs correlate with the unidentified 
IR emission features commonly attributed to polycyclic aro¬ 
matic hydrocarbon molecules. 

The Na° velocities within the LMC cluster around 
~ 280 km s“^, down to ~ 240 km s“^. The red “LMC” 
Ca"^ component is mostly found at velocities ~ 280 km s“^, 
i.e. associated with the bulk of the Na°. We can thus con¬ 
clude that this absorption arises within the LMC. The blue 
“LMC” Ca"*" component is mostly found around ~ 220 km 
s“^, i.e. outside the Na° range. In contrast, it is the blue 
“SMC” Ca"*" component which coincides with the ~ 140 km 
s“^ of the Na° within the SMC, whereas the red “SMC” 
Ca"*" component is found around 200-220 km s“^. 

The similarity between the blue “LMC” and red “SMC” 
Ca^ kinematics, and their weaker spatial resemblance to the 
Magellanic Clouds, suggests that this complex is not inter¬ 
nal to the Magellanic Clouds, but instead arises within the 
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Figure 9. Maps in the equivalent width of the 5780 and 5797 A DIBs, and their ratio, for the Galactic foreground in the direction of 
{top row:) the SMC and {bottom row:) the LMC. The locations of the target stars are indicated with the small white squares (detections) 
and magenta diamonds (non-detections). 


Milky Way Halo. De Boer, Koornneef & Savage (1980) sug¬ 
gested that UV absorption seen in front of the LMC at 220 
km s“^ has its origin within coronal gas associated with the 
LMC. Call absorption was already detected down to such 
low velocities towards 30 Doradus by Blades & Meaburn 
(1980), and a ~ 200 km s”*^ component was seen in Hi 
towards the LMC and Magellanic Stream by McGee, New¬ 
ton & Morton (1983). Intermediate- and high-velocity clouds 
such as those populating the Halo are generally transparent 
in Nal and show up much better in Call (Smoker et al. 2015) 
- hence perhaps why we noticed additional Ca ll components 
but not in Nal. The fact that we detect this kinematic com¬ 
ponent even towards the SMC suggests that the gas is not 
solely connected with an LMC corona but perhaps forms 
part of a pan-Magellanic halo, possibly associated with the 
putative common dark matter halo invoked by Bekki (2008). 


6 SUMMARY OF CONCLUSIONS 

We conducted a high signal-to-noise spectroscopic survey of 
666 early-type stars within the SMC and LMC, and mea¬ 
sured the absorption strength of the 5780 and 5797 A DIBs 
and of CallK and Nal. We examined correlations between 


these different ISM tracers and constructed for the first time 
maps covering much of these galaxies. Our findings are: 


• The 5780 A DIB is found in more strongly irradiated 
and/or hotter environments than the 5797 A DIB, but it 
too disappears when conditions become too extreme. Both 
DIBs are depleted in dense molecular clouds. 

• The 5780 and 5797 A DIBs have similar strength w.r.t. 
Nal in the SMC, LMC and Milky Way. This suggests that 
the abundance of their carriers is directly proportional to 
overall metallicity. 

• At low column densities, structure is seen at sub-pc 
scales in all the above tracers, and their correlations exhibit 
large scatter. Good correlations between different tracers 
seen in denser columns are probably due at least in part to 
the averaging effect of multiple, distinct structures along the 
sight-line even if different conditions prevail within them. 

• A common kinematic component around 200-220 km 
s“^ is seen in the CallK line towards both Magellanic 
Clouds, suggesting it arises in a pan-Magellanic halo. 


© 2015 RAS, MNRAS OPP.IHflSl 




































DIB maps across the Magellanic Clouds 13 


0.00 0.05 0.10 0.15 0.20 0.25 0.30 


equivalent width (Angstroem) 



1.2 1.1 1.0 0.9 0.8 0.7 

Right Ascension (hour) 


0.0 0.2 0.4 0.6 0.8 


equivalent width (Angstroem) 



No I D, 


1.2 1.1 1.0 0.9 0.8 0.7 

Right Ascension (hour) 



Figure 10. As figure 9 but for CaiiK and NaiD 2 . 
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APPENDIX A: SPECTRAL TYPES 

We performed a search in the literature for know spec¬ 
tral types, thereby mostly relying on the use of SIMBAD 
(Wenger et al. 2000); a summary of the pertinent sources for 
these spectral types is presented in table Al. These spectral 
types have been included in our catalogue, and comprise 174 
SMC and 61 LMC targets. Their distribution over spectral 
type and luminosity class is presented in figure Al. The more 
extensive coverage of the SMC includes also more later (A) 
type stars, whilst the known spectral types for LMC sources 
tend to include more O-type classes. Most of the targets with 
known spectral type are mid-B (SMC) or early-B (LMC). 

The question arises concerning the contribution from 
stellar photospheric absorption to the measurements of the 
calcium and sodium lines. Indeed, equivalent widths of the 
CallK line can reach 0.1-0.2 A at solar metallicity (Mihalas 
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show the expected scaling with metallicity - but this can 
also be explained if it were of purely interstellar origin. 

^ This all appears to be rather reassuring and lends cred- 

S ibility to the analysis of the ISM absorption. Spectral typing 
of the remainder of the target sample will be the subject of 
a forthcoming paper. 


Figure Al. Distribution over spectral type and luminosity class, 
for the (top:) SMC and (bottom:) LMC. 
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spectral type (02-A7) 

Figure A2. Distribution over spectral type and luminosity class, 
for the (top:) SMC and (bottom:) LMC. 


1973; Keenan et al. 1988). Figure A2, however, exhibits no 
dependence of the CallK absorption on spectral type, ex¬ 
cept perhaps for the A type stars that seem to have larger 
equivalent widths (but that are few in number). The equiv¬ 
alent width of the stronger (blue) kinematic component is 
mostly 0.2-0.4 A and thereby larger than what is expected 
for the stellar photosphere of (metal-poor) O- and B-type 
stars; the other (red) kinematic component is often much 
weaker but the same insensitivity to spectral type holds true. 
The same conclusion is reached for the NalD 2 (blue) and 
Di (red) lines (bottom panel in figure A2). 

We can also remark that the CallK absorption is of 
similar strength in the SMC and LMC, whereas we would 
have expected a clear difference if it mainly originated in 
the stellar photospheres due to the considerable difference 
in metallicity between the young populations of the SMC 
and LMC. The NalD absorption, on the other hand, does 


© 2015 RAS, MNRAS OOO.fTHT^ 






















